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We study the twist-3 fragmentation functions H and H, by applying a spectator model. In the 
calculation we consider the effect of the gluon rescattering at one loop level. We find that in this 
case the hard-vertex diagram, which gives zero contribution to the Collins function, does contribute 
to the fragmentation function H. The calculation shows that the twist-3 T-odd fragmentation 
functions are free of light-cone divergences. The parameters of the model are fitted from the known 
parametrization of the unpolarized fragmentation Di and the Collins function Hi . We find our 
result for the favored fragmentation function is consistent with the recent extraction on H and H 
from pp data. We also check numerically the equation of motion relation for H, H and find that 
relation holds fairly well in the spectator model. 
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I. INTRODUCTION 


The Collins effect [ljj has played an important role in the understanding of single spin asymmetries (SSAs) in various 
high energy processes, such as semi-inclusive deep inelastic scattering (SIDIS), hadron production in pp collision, and 
e + e _ annihilation into hadron pairs. The mechanism can be traced back to the so called Collins fragmentation 
function jj, denoted by -ffjS which is a transverse momentum dependent (TMD) nonpertubative object entering 
the factorized description of hard processes. It originates from the correlation between the transverse momentum 
of the fragmenting hadron and the transverse spin of the parent quark. Different from the ordinary unpolarized 
fragmentation function D\, the Collins function is time-reversal-odd and chiral-odd. The extraction of the Collins 
function has been performed in Ref. jgj , and in Ref. Q by considering TMD evolution. 

For quite some time it was believed that the dominant contribution to the transverse SSA for hadron production 
in pp collision comes from the the Qiu-Sterman function Tp(x,x) 3,1 , which can be related to the transverse- 

momentum dependent (TMD) Sivers parton density fi~{x,p\) jb;]: Tp(x,x) = — f cPptjS- fir( x ,P t)\sidis- The later 
one also contributes to the Sivers SSA in semi-inclusive deep inelastic scattering (SIDIS) under the TMD factorization. 
However, a recent study (7] showed that the function Tp{x,x) extracted from p^p —> hX does not match the sign 
of the Sivers function fitted from SIDIS data. This is the so called “sign-mismatch” puzzle. It was suggested Q 
that the twist-3 fragmentation contribution may be important for the SSA in pp collision, and could be used to solve 
the puzzle. This was further confirmed by a phenomenological analysis Q on SSA of inclusive pion production in 
pp collision [M3 within the collinear factorization, showing that the fragmentation contribution combined with 
the Tp(x,x) extracted from SIDIS data can well describe the SSAs in p^p —> ttX. In this framework, three twist-3 
fragmentation functions, H(z), H(z) and Hpjj(z, z±), participate. The first one corresponds to the first moment of 
the TMD Collins function and has been applied to interpret the SSA in pp collisions in previous studies BEl. 
The second one appears in subleading order of a 1/Q expansion of the quark-quark correlator, while its TMD version 
H(z, kp) is also a twist-3 function. The function Hpu(z, Z\) is the imaginary part of Hpp{z , Zi), which involves the 
F-type multiparton correlation 30,0. The three functions are not independent, as they are connected by the 
equation of motion relation 


dz 


H{z ) = -2 zH(z) + 2z 3 / -^PV-- -Hf v (z, Zl ) = -2 zH{z) + H(z). 


(1) 


In the last equation we have used H{z) to denote the “moment” of Hp rj (z, Z\). The function H might also contribute 
to the sin <j>g SSA in SIDIS through the coupling with the transversity distribution 0|. 
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Except for H, currently the quantitative knowledge about the other twist-3 fragmentation functions mainly relies 
on the parametrization in Ref. |9j. These fragmentation functions not only play crucial role in the understanding of 
the SSA in pp^ —> hX process, but also give significant contribution to the SSAs in single-inclusive leptoproduction 
of hadrons: IgX —> hX collision E3- The fra gme ntation contribution at the twist-3 level also enter the description of 
the longitudinal-transverse spin asymmetry [18|] in the process —» hX. Therefore, it is important to perform 

further theoretical and model study to provide information of H and H complementary to the phenomenological 
analysis. Besides, the function H{z) also encodes interesting information regarding the quark-gluon-quark correlation 
during the parton fragmentation. In this work we will study those fragmentation functions from the model aspect. 
Particularly, we will perform a calculation on the function H and H for the first time, using a spectator model. This 
model has been applied to calculate the Collins function for pions [l9l - l24| and kaons [24j], by considering the pion 
loop, or the gluon loop. In our calculation we will incorporate the effect of the gluon loop. We first calculate the 
TMD function H(z,k and H(z,k f.). The corresponding collinear functions are obtained by integrating over the 
transverse momentum. 


II. SPECTATOR MODEL CALCULATION OF H AND H 

Here we setup the notations adopted in our calculation. We use k and Ph to denote the momenta of the parent 
quark and the final hadron, respectively. We also apply the following kinematics: 

k = (k-,k + ,k T )=(k-,^-^,k T ), P h = (P^,P+,0 T )=\zk-,^ z ,0 T ,j, (2) 

where the light-front coordinates aX = a • rP have been used, kr denotes the momentum component of the quark 
transverse to the two light-like vectors n , and z = P^ jk~ is the momentum fraction of the hadron. The transverse 
momentum of the hadron with respect to the parent quark direction is given by Kt = — zkx■ 


A. Calculation of H up to one gluon loop 


The fragmentation function H(z 1 k^) can be obtained from the following trace 


^e^PH(z,kr) = ^Tr[A(z, k T )ia a P -/ 5 ], 

J h 1 

where A(A, fc-r) is the TMD correlation function that is defined as: 

1 




^yf 1 eIfe ' ? <°l ^(01^ ^[^(0) Z'(o+,00 + ) Z '(~^00 T ) |0> 


( 3 ) 


• ( 4 ) 


£-=o 


Here U? ^ denotes the Wilson line running from a to b at the fixed position c, to ensure the gauge invariance of the 
operator. In the spectator model, the tree level diagrams lead to a vanishing result because of lack of the imaginary 
phase. To obtain a nonzero result one has to go to the loop diagrams. In one-loop level there are four different 
diagrams (and their hermitian conjugates) that may contribute to the correlator A(z,kj,), as shown in Fig. [I] These 
include the self-energy diagram (Fig. [lji), the vertex diagram (Fig. [!]:>), the hard vertex diagram (Fig, [lj:), and the 
box diagram (Fig. [lji). They have also been applied to calculate the Collins function in Refs. [23l |24|. 

We will focus on the the favored fragmentation function, i.e. the fragmentation of u —> 7r + . In this case the 
expressions for each diagram in Fig. [T| are as follows: 


A( n )(z, kr) 


A(i,)(2, k T ) 


= i 


= % 


4 C F a s (ty + m) ! u 

2(27t) 2 (1 - z)P~ (fc 2 - to 2 ) 3 9qhlZ 1 * 
f d 4 l Y ( - / + m) ( # + m) 

(27t) 4 {(k — l) 2 —m 2 + ie)(l 2 + ie) ’ 

4 C F a s (1f,+m) (1L 

2(27t) 2 (1 - z)P^ ( k 2 - m 2 ) 2 9qhlb 1 * 
f dH 7 /J (#- tfh ~ / + m s )g qh j 5 ( If. - 
(27t) 4 ((fc - Ph — l) 2 — Tn 2 + ie)((k - l) 


tfh + m s )g q h"lb{ H + m) 


ifh + m s ) 

- / + to)7 m ((If + m) 

2 — to 2 + ie)(l 2 + ie) ’ 


( 5 ) 


( 6 ) 
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FIG. 1: One loop level diagrams utilized to calculate the correlator in the spectator model. The double lines in (c) and (d) 
represent the eikonal lines. The hermitian conjugations of these diagrams, which we have not shown here, also contribute. 


\c){z,k T ) = i 


A(d)(z, k T ) = i 


4:CfOI S 


2(2tt) 2 (1 -z)P- k 2 ~m 
f d A l 7 + (\f — / + to) 


(U + to) , „ . , „ , 

- 2 9qhlb {$- Ifh + m s )g qh 75 ( # + to) 


(27t) 4 ((fc — l ) 2 — to 2 + ie){—l~ ± ie){l 2 + ie) ’ 
4:C F a a ( ^ + to ) 

- 5 - {If - Ifh + m s ) 


2 ( 2 tt ) 2 (1 - z)P* k 2 - ' 


d A l 

(27r) 4 ((fc — Ph — l) 2 — to 2 + ie)((k — l) 2 — to 2 + ie)(-l~ ± ie){l 2 + ie ) 


7+( - If h - / + TOg) g 9 /i7 5 ( ty - / + to) 


( 7 ) 


( 8 ) 


Here g q h is the coupling of the quark-hadron vertex, to the mass of the quark in the initial state, and m s the mass of 
the spectator quark. In Eqs. 0 and (H|) we have applied the Feynman rules for the eikonal lines. 

In the calculation of T-odd functions, one should utilize the Cutkosky cut rules to put certain internal lines on 
the mass shell to obtain the necessary imaginary phase. For T-odd fragmentation functions, only the cuts through 
the gluon line and the intermediate quark line inside the loop give rise to the result. This corresponds the following 
replacements 


1 

l 2 + ie 


—> —2ni5(l 2 ), 


1 

(,k — l) 2 + ie 


—> —27 HS((k — l) 2 ). 


( 9 ) 


Here the cuts thr oug h the eikonal lines do not contribute. This directly links to the universality of the TMD fragmen¬ 
tation functions |25H28l ] . which has been verified intensively in literature. Another issue that should be addressed is 
the choice of the quark-hadron coupling g q h. When choosing the point-like coupling, there is a divergence appearing 
at large k F region in the calculation of the collinear fragmentation function: 


H{z) 


d 2 K T H{z , fc|) = J d 2 k T H(z, *4). 


( 10 ) 


In the literature two different approaches have been applied to regularize this divergence. One strategy is to adopt a 
cut on k F by putting an upper limit fc™ aa: , The other is to choose a form factor for g q h which depends on the quark 
momentum. Here we will utilize the second approach. Follow the choice in Ref. [24| , we adopt a Gaussian form factor 
for the coupling, 


_ fc 2 

e a^" 

dqh^gqh - ( 11 ) 

Z 

where A 2 has the general form A 2 = A 2 /(z“( 1 — z)^). The A, a, and /? are the parameters of the form factor that will 
be determined in the next section. The advantage of the choice in Eq. [TT] is that it can also reasonably reproduce [2dj 
the unpolarized fragmentation function. 

In Eqs. (HI) or (|5|). in principle one of the form factors should depend on the loop momentum l. Here we will drop 
this dependence and merely use k 2 instead of (k — l) 2 in that form factor to simplify the integration. The same 
choice has also been adopt to calculate the Collins function [24}, which is a leading-twist fragmentation function. For 
the subleading-twist T-odd functions the situation is more involved. As shown in Ref. [ 29 }, the calculation of T-odd 
twist-3 TMD distributions suffers from a light-cone divergence. In phenomenological studies the divergence has to be 
regularized |29l l30[ by introducing form factors, explicitly depending on loop momentum. However, as we will show 
later, we find that in the case of twist-3 fragmentation functions, the calculation is free of this light-cone divergence. 
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The reason behind this distinction is that the kinematical configuration contributing to T-odd fragmentation functions 
is different from that to the T-odd distribution functions. 

After performing the integration over l using the cuts in Eq. O we organize the expression for H(z,k^) as follows 


H(z, fcy) — 


2 a a g*C F 


1 


(27r) 4 z 2 (l — z) Mh(k 2 — m 2 ) 

The four terms in the bracket of the right hand side of (fT2l) have the forms 


(H {a) (z,k^) + H {b) (z,k^) + H {c) (z,k^) + H [d) (z,k^)) 


H(a) (A ^t) — 


_^_( 3 _ EL) (J fc» - m 2 + (1 - 2/*K)J 1 


H( b ) (z, k%) = ( k mfl + ™ s h - m s / 2 ) (k 2 - m 2 s + (1 - 2 /z)m 2 h ), 

V \{m h ,m s ) ) 

H( c )(z, k^) = —((it is — m)(k 2 — mm s ) + mmh)Ii/(k 2 — m 2 ) — (m s — m + zm)l 3 ,k~ , 

H( d) (z, k?) = -m + zm)(\(m s ,m h ) + ((1 - 2 z)k 2 +m 2 h - m 2 s ) ( k 2 -m 2 s + ( 1 - 2 /z)m 2 h ) 

— zm (k 2 — m 2 + (1 — 2 /z)m 2 1 )) / 2 — I 2 ((m s — m)(k 2 — mm s ) + mm 2 ) + (m s — m + zm)l 3 k~ 
The functions represent the results of the following integrals 

Ii = j d 4 M(Z 2 )5((fc — l) 2 — m 2 ) = ^2 (^’ 2 — w 2 ) > 


h = / cH 


S(l 2 )6((k-l) 2 -m 2 ) 


(k-P h - l) 2 - 


2A (mh,rn s ) 


In 1 + 


2y/\(m h ,m s ) 


k 2 — m 2 h + m 2 + yj\(m h ,m s ) 


h= f dil 5(l 2 )5((k-l) 2 -m 2 ) 

J —l is 


( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 
(19) 


with A (to/!, m s ) = (k 2 - (mu + m s ) 2 )(k 2 - (run - m s ) 2 ). 

We would like to point out that the quark-photon hard-vertex diagram gives nonzero contribution to H(z, kr), as 
shown in Eq. [T5l This is different from the calculation of the Collins function H j 1 , in which case the contribution 
from the hard-vertex diagram vanishes (24|. We note that this is because the Dirac structure of H(z,k f.) appearing 
in the decomposition of the correlation function A (z, hr) is different from that of the Collins function. The sum of 
H( c )(z,kr) and H^(z,kr) can be cast into 


H( c+ d)(z,kl) = -m + zm)(X(m s ,m h ) + ((1 - 2 z)k 2 + m 2 h - m 2 s ) (k 2 - m 2 + (1 - 2 /z)m 2 h ) 

- zm (k 2 - m 2 + (1 - 2 /z)m\)) h - + h'j ((m s - m)(k 2 - mm s ) + mm 2 ) , 


( 20 ) 


where the terms containing I 3 cancel out. As we can see, the final result of H(z, kj,) in Eq. (fTSl) is free of the light-cone 
divergence. 


B. Calculation of H with gluon rescattering 


The fragmentation function H(z,k^) originates from the quark-gluon-quark (qgq) correlation 0 , 111 : 

/>*,n, f i+ 

(2 tt )3 J 


A ^(zM) = 


sh 1 


2 zN c 


b 

J x 


/ ±00+ 


dl l +U f^+,v+) 


x gF- a ( V )U^ + ^(0\Ph;X)(P h ; X\^(0)U^ + ^ +) U$ T ^ t) |0) 


? 7+=£+=0 


( 21 ) 
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FIG. 2: Diagram relevant to the calculation of the qgq correlator in the spectator model 


where F M1/ is the antisymmetric field strength tensor of the gluon. Using the identity 


r£ + r °° + 

/ dg + =± dg + 9(±£ + =F V + ) 

J ±00+ OO + 


■s 


drf 


1 1 \ 


d I -- 

z Zl 


(i-i) 


( 22 ) 




with 6 is the Heaviside function, we can rewrite the qgq correlator as 


A A(z,k T ) = 


sh 1 




2zN r 


d^ + d 2 ^rdg^ 

(27T) 4 


di 1 —^ 6 

Z Zl 


x(0\igFZ a (v)m\Ph;X)(P h ;X\m\0) 


- —-— ie 

Z Z l 




ri+=£+= 0 
TjT—^T 


(23) 


Here we have suppressed the Wilson lines for brevity. In Eqs. (1^1) and (E3l) we use \/z — \/z\ to denote the momentum 
fraction (along the minus light-cone direction) of the gluon with respect to the final state hadron, following the 
notations Ref. @ . Thus 1/ Z\ gives the momentum fraction of the quark correlated with the gluon. 

The fragmentation function FI can be extracted from the correlator A“(z, hp) by the following projection: 


; Tr[A°^(z, k T )cr a } = H(z, k?) + iE(z, kj.) ■ 


(24) 


The integrated fragmentation function H(z) = z 2 f d 2 krH(z,kp) is related to the collinear twist-3 fragmentation 
function F[p U (z, z\) by 


r OO 7 -1 

H(z) = 2z 3 / Zl ). 


(25) 


where H% r ( z, z 1 ) is the imaginary part of Hpu(z, Z\) that appears in the decomposition of the F-type collinear 
correlator |gj, [lj| 


J x 


d t + 


h 


27T / 27T 


(0\igF I a (r 1 + )^+)\P h -,X)(P h ;X\m\0) 


= M h 


-°L a /3 + 75 H FU {Z, Zl) 


(26) 


The diagram used to calculate the fragmentation function H in the spectator model is shown in Fig. [21 which 
represents a qgq correlation. The left hand side of Fig. [2] corresponds to the quark-hadron vertex (P^; X|t/>( 0)|0), 
which has the following form in the spectator model 


U{Px){ils) 


i(lf + m) 
k 2 — m 2 ’ 


(27) 














6 


with Px denoting the momentum of the spectator quark. The right hand side of Fig. [2] corresponds to the vertex 
(0\igF± a (g + )il)(£, + )\Ph', X), whose expression can be given in a similar way. The differences are that one should 
consider the field strength tensor F a/3 , as denoted by the circle at the end of the gluon line in Fig. [2] Its Feynman 
rule (on the right hand side of the cut) is given by i(q a g^ p — q^g ap )5 a b, with p and b the indices of the gluon line. 
Thus, we can write down the expression for the qgq correlator as: 


K a A { Z ,k T ) 


4 Cfol s 1 

*2(2-7t) 2 (1 — z)P^ k 2 — to 2 

d 4 l (/~g^ - It 9~ p ){ # ~ / + m ) g q hlb {If - tfh ~ / + rn s ) 7m ( - rfh + m s ) g qh y 5 ( # + m) 

(2-7r) 4 (— 1~ ± is)((k — l) 2 — m 2 — is)((k — Ph — l) 2 — m 2 — is)(l 2 — is ) 


where we have used the replacement 


i-i| p. 

Z Z\ 


(29) 


According to Eqs. (12411 and (|28|l. the contribution to H comes from the imaginary part of sub-diagram shown on 
the right hand side of the cut in Fig. [2] In order to do this, again one needs to apply the Cutkosky cutting rules 
to integrate over the internal momentum l, that is, to consider all the possible cuts on the propogators appearing in 
Eq. (1281) . However, only the cuts on the gluon line and the fragmenting quark survive, as shown by the short bars 
in Fig. [2] Other combinations of cuts are kinematically forbidden or cancel out each other. In particular, the total 
contribution from the pole of the eikonal propagator is zero. To demonstrate this, we consider two different cases. 
The first case is to take the poles of l/(— 1~ ± is) and 1 /(l 2 — is), therefore, It has to be zero. This yields vanishing 
H since there is a factor l~ g^jA ~ ^t9~ p 411 the numerator of Eq. (1551) . The second case is that one applies the cut on 
l/(—l~ ± is) and 1 /((k — l) 2 — m 2 — is), or on l/(— 1~ ± is) and 1 /{(k — Ph — l) 2 — ml — is). However, these two 
contributions cancel out each other. This is because the pole positions for l + from the propogators 1 /(( k—l) 2 —m 2 —is) 
and l/((— 1~ ± is) and 1 /((fc — Ph — l) 2 — ml — is) are on the same half plane, which means that the integration over 
l + vanishes with the delta function S(l~) (since k~ — P^ > 0) 


dl+ 


1 


27t ((fc — l) 2 — m 2 — is){(k — Ph — l) 2 — m 2 — is) 
f dl+ 1 


2tt (2jfe-(jfc+ -l+) + -*e)(2(Jfe- - P^){k+ - P+ - 1+) +- is) 


= 0 


(30) 


Therefore, we will again apply the cutting rules given in Eq. ([21) to perform the integration over l, and the factor 
1/(1 /z — 1/zi ± is) will take the principal value, as also shown in Refs. [32, 33|. This means that H is process 
independent in the spectator model, similar to the Collins function and H. The final result for H has the form 


H(z, k ^) — 


Oi s g, 


qn 


(27T) 4 


C F 


e a-* 


1 


2 2 (1 -z) M h {k 2 -m 2 ) 


— -4(m s — m)fcy 


+ [m s — m + zm) [A(k 2 + k^) + BM^/z — I\/z — ( k 2 — m 2 )^/-] 
+ [( k 2 - mm s )(m s — m) + mm 2 ] [ h/(2zk 2 ) + A/z + B\ > . 


Here A and B denote the following functions 

h 


A = 


B = - 


which appears in the integration 


A (m h ,m s ) 
2 k 2 


( 2 k 2 ( k 2 -m 2 s - m 2 h ) + (k 2 + m 2 h - m 2 s )^j , 


A (m h ,m t 


. 7l (l + k2 + m 2 s -ml h 


4l l p 8(l 2 ) S((k — l) 2 — m 2 ) 
[k- P h - l) 2 - m 2 


dr l 


= Ak p + BP p . 


( 31 ) 

(32) 

(33) 

(34) 
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m s (GeV) 

A (GeV) 

gqiT 

m (GeV) 

a 

0 

0.53 

2.18 

5.09 

0.3 (fixed) 

0.5 (fixed) 

0 (fixed) 


TABLE I: Fitted values of the parameters in the spectator model. The values of the last three parameters are fixed in the fit. 


III. NUMERICAL RESULT 

In this section we present the numerical result for the fragmentation functions H and H. To this end the values of 
the parameters in the model have to be specified. In Ref. [2J] the parameters of the model were determined by fitting 
the model result of unpolarized fragmentation function D\(z) with the Krezter parameterization j[34| of D\(z). The 
parameters were then used to make prediction on the Collins function. In this paper we will obtain the parameters by 
fitting simultaneously the model calculations of the unpolarized fragmentation function and the Collins function with 
the known parameterizations of them, since the Collins functions have been extracted and are well constrained by the 
e + e _ annihilation data and the SIDIS data. Specifically, we will use the half-fc-p moment of the Collins function 

Ht {1/2 \z) = z 2 f d 2 k T ^Ht{z,k 2 T ) (35) 

J 2 m h 

in the fit. 

For the theoretical expressions of D\ and H ^~, we use the calculation in the same model, which has already been done 
in Ref. [24[ 1 . For the parameterization of D \, we will adopt the DSS leading order set (35| . For the parameterization 
of the Collins function, we apply the recent extraction By Anselmino et.al. [§]. We note that in Ref. @|, the DSS 
fragmentation function is also used to extract the Collins function. 

Our model calculation is valid at the hadronic scale which is rather low, while the standard parametrization of D\ 
is usually given at Q 2 > 1 GeV 2 . Therefore we extrapolate the DSS D\ fragmentation to that at the model scale 
Q 2 = 0.4 GeV 2 in order to perform the fit. For the same reason, the Collins function should be evolved at that scale 
for comparison. However, the evolution of the Collins function is rather complicated Q HH, ;33]- In the extraction of 
the Collins function in Ref. @|, the authors used the assumption that the Collins function evolves in the same way of 
Di(z). The same assumption has also used in Ref. [241] For consistency we will use this assumption since in the fit we 
use the parametrization of Collins function from Ref. |2| . 

In Table. H] we list the fitted values of the parameters in the model. In the left panal of Fig. [3l the curve (the 
solid line) vs z for the unpolarized fragmentation function Di(z) at the model scale Q 2 = 0.4GeV 2 is compared with 




FIG. 3: Unpolarized fragmentation function D\(z) (left panel) and the half moment of the Collins function (right panel) vs 
z for the fragmentation u —f 7r + at the model scale Q 2 = 0.4GeV 2 . The parameters are fitted to the parameterizations in 
Refs. [35} and Q. The result in Ref. [24] (dashed lines) is also shown for comparison. 


1 We recalculate the Collins function and find that our result does not exactly agree with the result in Ref. Vt |. For completeness we 
present our result for II ] in the Appendix. 








































FIG. 4: Left panel: The twist-3 fragmentation functions H(z) and H(z) vs z, plotted by the solid line and the dashed line, 
respectively. Right panel: H(z) compared with —2 zH(z) + H{z) in the spectator model. 


the curve (dotted line) from the DSS parameterization. We also show the result (dashed line) calculated from the 
parameters fitted in Ref. [24j . In the right panal of Fig. 0 we display the fitted curve for h[ 1 ^ 2 \z) and compare it 
with the parametrization of Ref. Q. 

In the left panel of Fig. [I] we plot our prediction on H(z) and H(z ) using the parameters in Table. [Q We present 
the result at the model scale Q 2 = 0.4GeV 2 , We find that the sign of the favored H(z ) is negative and its magnitude 
is sizable. This is consistent with the extraction in Ref. Q, where a negative H(z) for the favored fragmentation is 
given. For the function H(z), we find that the result is nonzero and has a minus sign, in Ref. ji|, a similar result is 
also hinted by the fit on Hpjj(z, zi), which contribute substantially to H(z ) through Eq. [L] 

According to Eq. [I] the three twist-3 fragmentation function should satisfy the equation of motion relation, which 
is a model independent result derived from QCD. However, From Eas. dT^l) . (1311) and (1M1) . one can not find out an 
obvious relation among them since in the spectator model they are calculated from different diagrams. Thus we 
numerically check the relation m and show the the comparison between H(z) (solid line) and —2 zH(z) + H(z) 
(dashed-dotted line) on the right panel of Fig. [4] We find that the two curves are close, which indicates that the 
relation holds approximately in the model, therefore it provide a crosscheck on the validity of our calculation. 


IV. CONCLUSION 

In this work, we studied the twist-3 fragmentation function for H and H in a spectator model. We first calculated 
the TMD functions H(z, kp) and H(z , fcf.), and then we obtained the corresponding collinear functions by integrating 
over the transverse momentum. In our study we considered the gluon rescattering effect and found that the hard- 
vertex diagram gives nonzero contribution to H. Using the parameters fitted to the known parameterizations of D\ 
and Hi simultaneously, we presented numerical results of H and H. We found that our results agree with the recent 
extraction from the SSA in pp collision. We also tested the equation of motion relation among H(z), H(z ) and H(z), 
the numeric result shows that the relation approximately holds in our calculation. Our study may provide useful 
information on the twist-3 fragmentation function complementary to phenomenological analysis. 
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Appendix A: Results of the Collins function 

Here we present the model result of the Collins function (24j 


H\{z, fcy) — 


2a a gLC F 


M h 


(27r) 4 z 2 ( 1 — z) ( k 2 — m 2 ) 


(^i (a) (a &t) + H\(b) ( z > ^t) + Hi(d) i z i ^r)) (36) 


The three terms in the brackets correspond to the results from Fig. [TJi, Fig. [[):>, and Fig. QJ1, respectively. In our 
calculation we find that those terms have the form 


H 1 (a)(-> fc T) = 


TO 


( k 2 — TO 2 ) 


3 - 


TO 

fc 2 ' 


H\( b ){Zi fcr) = 2 m s I 2 - 2(m s - to) 


h 

m 2 — to? — k 2 


X(m h ,'i 


-h - 


Ak 2 m 2 
\(m h ,m s )Tr 


hk 


H^ d ,(z, k t) = 2 {-/ 34 (2zto + 2 m s - 2to) + / 2 [2zto (A: 2 - to 2 + M 2 (l - 2/z)) 

Z Z K/rj-i 


+ 2(to s — to) ((2z — l)/c 2 — M% + to 2 — zm(m + m s ))]}■ 
Here / 3 4 is the combination of two integrals 


hi = k (J 3 + (1 - z)(k 2 - m 2 )h) = 7rln 


vW^) 


TO S 


(37) 

(38) 

(39) 

(40) 


with 


/4 


d 4 l 


S(l 2 )S((k — l) 2 — to 2 ) 

(—h + ie){k — p — l) 2 — m 2 


(41) 


We find that in (1551) there is a new term proportional to m s — m that was not contained in Eq. (29) of Ref. 24 [. 
Also in Eq. (l39l) the coefficient of certain terms containing m s — to has a factor of 2 compared to Eq. (30) of Ref. 241. 
But our calculation returns to the results in Ref. [23} in the case m s = m and by setting the form factor to 1. 
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